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ABSTRACT 

We  hope  to  answer  one  of  the  most  fundamental  and  important  unsolved  ques¬ 
tions  in  chemistry:  how,  from  a  molecular  perspective,  do  diemical  reactions  in  solu¬ 
tion  actually  occur.  The  key  to  solving  this  long-standing  problem  is  to  understand 
the  molecular  dynamics,  i.e.  the  motions  of  the  atoms  and  the  forces  that  drive  them. 
We  have  already  developed  theoretical  techniques  and  computational  procedures 
involving  specialized  computer  hardware  needed  to  calculate  the  molecular  dynamics 
for  many  chemical  reactions  in  soiutirm.  From  the  dynamics  we  have  derived  the 
interface  for  experimental  verification,  namely  transient  electronic,  infrared,  and 
Raman  spectra  as  well  as  x-ray  diffractimt,  all  of  which  are  potentially  observable  man- 
ifestetions  of  the  atomic  motions  during  the  reaction.  We  have  tested  our  approach 
on  the  simple  inorganic  I2  photodissociation  and  solvent  caging  reaction.  The  agm- 
ment  between  molecular  dyruunics  based  theory  and  expnimental  picosecond  transient 
electronic  absorption  spectrum  as  a  function  of  solvmt,  time,  and  wavelength  is 
sulUently  dose  as  to  indicate  that  for  the  first  time  we  are  discovering  at  least  part  of 
the  moleciilar  dynamics  by  which  a  real  solution  dmnical  reaction  takes  place. 
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To  be  puMWied  in  Prooeedims  oT  the  Imemeiionel  Conference  on  Photothemistry  and  PhoioMolosy  (Aten- 
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I.  INTRODUCTION 

The  vast  majority  of  the  reactions  of  interest  to  chemists  occur  in  sdution,  including  most  of 
inorganic,  organic,  industrial,  and  bio-chemistry.  It  is  thus  ironic  that  there  is  not  yet  a  single  solu¬ 
tion  chemical  reaction  is  understood  from  the  microscopic  point  of  view  of  the  atomic 

motions,  i.e.  the  molecidar  dynamics,  by  which  it  occurs.  Therefore,  one  of  the  fundamental  and 
important  questions  of  chemistry  remains  unsolved.  Answering  this  question  is  important  not  just 
from  the  viewpoint  of  pure  research,  but  also  because  an  understanding  of  how  solution  reactions 
actually  occur  can  be  expected  to  lead  to  an  increased  ability  to  control  such  reactions  for  useful 
ends. 

We  believe  that  there  is  a  general  approach  which  can  lead  to  the  discovery  of  the  molecular 
dynamics  for  a  wide  variety  of  inorganic,  orguiic,  and  biochemical  solution  reactions.  Such  a  gen¬ 
eral  approach  requires  both  theory  and  experiment,  as  only  theory  can  provide  a  complete  view  of 
the  molecular  dynamics  and  only  experiment  can  sufficiently  test  the  theoretical  predictions  that  we 
should  have  confidence  in  them.  Our  approach  has  three  parts.  First,  we  have  developed  the  com¬ 
putational  techniques  to  calculate  the  molecular  dynamics  for  reactant-solvent  systems  involving 
hundreds  to  thousands  of  atoms.  This  has  been  accomplished  by  the  use  of  specialized  computer 
hardware,  an  array  processor,*' 2  allowing  us  to  run  approxinuitely  35  times  faster  than  on  the  DEC 
VAX  780  we  also  use.  Second,  we  need  an  interface  between  theory  and  experiment,  a  set  of 
phenomena  which  closely  reflect  the  atomic  motions  by  which  reactions  occur  and  whidi  equally  can 
be  computed  theoretically  and  measured  experimentally.  Our  choice  for  this  interface  is  transient 
spectra.  Therefore  we  have  developed  the  theoretical  and  computational  techniques  to  calculate 
infrared,’’^  Runan,^  and  electronic^’*  spectra,  as  well  as  x-ray  diffraction*  from  molecular  dynam¬ 
ics.  Third,  we  have  begun  to  experimentally  measure  transient  spectra  from  chemicai  reactions  on 
the  picosecond  time  scale.  We  present  below  transient  electronic  absorption  from  a  diemical 
reaction,^*  as  both  theoretically  computed  from  molecular  dynamics  and  experimentally  measured. 

We  have  demonstrated  this  approach  for  the  chemical  reaction  of  I2  photocfissodation  in  vari¬ 
ous  solvents,  followed  by  solvent  caging,  atomic  recombinatitm,  and  the  decay  of  the  excess  energy 
in  the  newly  formed  vibrationally  hot  Ij  back  to  the  solvent.  We  have  computed  the  molecular 
dynamics  for  this  reaction  sequence  in  several  solvents,  and  from  the  dynamics  computed  the  time 
dependent  electronic  absorption^* *•**  and  Raman  spectra*'*  as  well  as  the  traroient  x-ray 
diffraction.*  To  check  the  dynamics,  we  have  constructed  a  picosecond  transient  electronic  spectrom¬ 
eter**' *2  and  measured  the  transient  etoctronic  absorption  for  this  reaction  in  a  variety  of 


solveatt.^**  Ite  nwamn  between  theofy  and  experiment  ai  a  ftmciioa  of  natura  of  aoivnt,  of 
time  aale,  and  of  ptobe  photon  enaiiy  ia  moat  iradfyiiit,  liviai  na  aome  maaon  to  hope  diat  for 
the  flnt  time  we  are  actualhr  diaooverinf  at  leapt  part  of  the  moleciilar  thmandre  by  whldi  a  real 
reaction  in  aolation  takaa  plaoe.  A  Iqrpotl^  b>  endphi  dip  wi^  vartatton  in  i»dppi|ipr  ^hoamk 
(as  well  as-apedfil-liiid^'amll  wiicir ye"' etafciW^bdBi ^dwowmbal^ lifa ■  oxpi^tteiiBy  among 
dUSBreot  solvents  is  oflbred  in  tenns  of  a  IMiP  ciaienae  picture  of  die  solvent  as  driven  by  the 
motions  of  the  reacting  molecules.  This  hypothesis  is  explored*<*  by  oomputtag  the  vdodty  spectra 
of  various  solvents  aadi  dhowIng^thBt  thay  do  indeed  soarelata  well  in  the  fnguenoy  range  of  interest 
with  the  observed  time  aesles  for  reaction  in  dUlHeat  aelveots. 


II.  M(>LECULAR  DYNAMICS 

The  computation  of  the  molecular  dyromics  of  solution  reactions  is  in  principle  easy.  We 
know  from  extensive  gas  phase  experience  that  dasskal  mechanics  will  usually  be  ad^uate.  Given 
the  potential  energy,  or  the  foroes,  linking  the  N  teoms,  and  a  set  of  initial  positions  r'^CO)  and 
momenta  p'*(0),  we  can  compute  the  traiectories  r’*(t),  as  shown  in  Fig.  1. 


Fig.  1.  TheofeHcd  udmique 
to  confute  abaorftkm  and 
scattering  spectra  fiom  moke- 
ular  dynamics.  The  intfec- 
lories  of  the  atoms  are  adcu- 
ktte^'^  from  classical 
mechanics  and  then  used  to 
confme  die  time  histories  of 
the  Interactien  trith  the  nalte* 
don  fieU,  Le.  die  dipole  mo- 
mem  fsU)  far  akseepdon  or 
the  pokolssdddiy  P(r)  frtr  Ra¬ 
man  scattering.  LiwiBr 
tesponse  iheoty,  estsembk 
amag/ag.  mtd  simple  quan¬ 
tum  conecdons  then  gjtee  the 
lafraitd.y*  or 

Ramatfr*  qpectra. 


In  practice  there  are  several  hundes  to  be  overeome.  One  is  die  aynthesia  of  the  best  poariMa 
potential  energy  surface  from  the  avaUabla  diaoretical  and  oxperimantal  inforiiMtion.  Wa  expect  in 
the  future  to  have  to  turn  increasing  for  gfdthmtea  to  oh  MHe  quantum  calculations,  a  direction  in 


wtrich  we  iMve  abea^r  tieiim  to  mova.*^  A  aacood  dMculty  to  tlia  ahaar  mmbar  of  mwiefical 
opentioM  wbidi  nniat  be  cufM  out  to  ooR^itte  a  ataiiatkaiy  nbaMe  deacfiftion  of  the  tj^Bandca, 
iq>  to  10*^  aritlmatic  opendiona  for  KT”  aeoooda  of  tbrnamica.  Iliia  wa  aolva  by  oae  of  an  tmv 
pmomotM  Other  dtaWentae  induda  the  proper  Aandliiit  of  hpiindaqr  coadhiena,#aubi<id  in  tea- 
aoaaHy  food  order  for  aochatfad  apadea,  but  not  yM  adiafacti^  aolved*^  for  lo^-finge  ICouiom- 
bic  interactiooa,  audi  aa  thoaa  involving  iona.  A 

While  we  do  not  want  to  minimiee  the  dibit  involved  in  oonqndnf  the  molecdtf  of 

aolution  leadkaia,  wa  beiiave  that  the  iargw  chaBenfe  now  liea  in  a|iiB^tly  teating  the  dyvamica 
againat  exparinad  ao  ttat  the  dynamics  can  dther  be  trusted  If  theory  |ipd  expoiment  agree,  or 
improved  if  thqf  do  not 


UI.  SPECTIA  nOM  MOLECULAR  DYNAMICS 

In  this  section  wa  brieQy  describe  some  of  the  types  of  spectra  which  we  are  able  to  compute 
from  solution  reaction  moleoular  dynamics,  and  whidi  we  wbh  to  experimentally  measure. 

Following  the  woric  of  Roy  Gordon  we  are  accustomed  to  dedudng  molecular  dynamics  fnm 
spectra.  Here  we  take  the  inverse  approadt  and  compute  spectra  from  molecular  dynamics. 
SpedflcaBy,  from  atomic  traiectories  we  conqwte  infrared,  nooresonance  Raman,  and  dectronic 
spectra,  as  wdl  as  x-ray  diffinKtion.  Fig.  1  illustrates  the  general  tedutique.  From  an  initial  set  d 
atomic  positions  and  vdodties  «d  the  forces  among  the  atoms  we  cdculate  from  Newton’s  second 
law  the  atomic  trdectories  r’*(t).  For  infrared  or  electronic  absorption  spectra  the  molecular  interac¬ 
tion  with  the  radiation  Add  is  through  the  dipole  moment  ft,  and  for  Raman  spectra  throu^  the 
pdarixability  tensor  1'  we  conqiute  ft(t)  and  P(t)  from  ft^^)  and  PCr**)  as  fundions  of  nuclear 
positions,  combined  with  the  traiectories  r'^ft)  as  functions  of  time.  The  technique  can  be  viewed  as 
an  expression  of  daasicd  linear  response  thewy  or  as  an  exerdse  in  classical  electromagnetic  theory. 
We  use  a  power  spectral  approach,  symbolized  as  DO,  talcing  advantage  of  fast  Fourier  analysis  with 
windowing  and  associated  corrections. An  average,  symbolized  by  <>,  of  the  spectrum  is  taken 
over  the  ensemble  appropriate  to  the  experimental  conditions  (for  example  a  constant  temperature 
and  density)  and  simi^  quantum  corrections  to  the  spectra  are  applied  where  needed. 

Iif/raitd 

In  order  to  compute  infrared  spectra,  we  cmnpute  the  molecular  dynamics,  i.e.  the  trajectories 
of  the  atcmis  r^(t),  and  substitute  the  trajectories  into  the  dipole  moment  function  to  produce 
the  time-varying  dipole  moment  function  from  which,  as  shown  in  Fig.  1,  the  infnured  spec¬ 
trum  is  computed. 

Fig.  2  shows  an  example  of  such  a  computed  infrared  spectrum,  for  the  gas  phase  Aindamental 
vibntiomd-rotational  band  contour  of  CX).  Note  that  alnoost  perfect  agreement  with  both  quantum 
medUHiical  theory  uid  experimental  measurement'^  is  achieved  by  an  essentially  dassical  mechani¬ 
cal  molecular  dynanrics  odculation  with  no  ac^ustable  parameters.  Elsewhere^  we  have  also  com¬ 
puted  pure  rotational  absorption  band  contours  and  have  extended  vibrational  and  rotational  infrared 
spectral  calculations  to  the  liquid  phase. 

B.  Raman 

To  compute  nonresonance  Raman  ^rectra  from  moleoular  dynamics,  we  compute  the  time 
vaiyiiii  pol^zability  matrix  P(t)  from  the  polarizability  matrix  P(r'*)  as  a  function  of  nucdear  posi¬ 
tions  and  the  trqject^es  r'*(t).  Fig.  3  shows  an  example  of  an  equilibrium  liquid  sute  calculation  for 
Nj  in  liquid  Argon  at  two  different  temperatures.’- Again  note  the  dose  agreement  of  molecular 
(hmandcs  calculation  with  experimoit,  without  the  need  of  adjusuMe  parameters. 


111.  2.  Eguttbrium  hifiwed 
apecmm  CO  gas  phae  fim- 
ibrnaital  vStrattonal- 

mtational  band  contour. 
Moleadar  tfynmrdcs  cdada- 
don  as  Matties?  accurate 
quantum  caktdatkm  as 
stni(fAr  Une,^  and  experimen- 
Ud  spectrtm^^  as  dotted  line. 
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Fl|.  3.  Eq^dUbtium  Raman 
spectra.  FUndamaUal  isotropk 
and  anlsoiroidc  stbradmud- 
rotadomd  Roman  spectra  for 
N2  In  Uqtdd  Ar,  at  two  tem¬ 
peratures.  Moleadar  t^m- 
ks  spectra  shown  as  solid 
llnef  dotsfiom  errperbrmt.^^ 
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Ekemtic 

By  wing  tlw  tmiritiM  dtoole  moraeot  faMtead  of  the  iMud  iBpote  OMNwatt  m 
cap  apto  uae  the  ttaiactoriaa  f"(t)  to  cikruhoe  the  thoe  variahon  of  |l•(t),  and  ftom  diia  die  elec* 
traaic  abooiption  apectnim^*'  in  a  panHel  wiv  to  the  teftarad  tdMoiy^  apacttoin.  Her  auflkien^ 
low  fiequency  motioo,  as  in  no  quantum  conecdon  is  neoesssiy.  4  shows  such  a  cowyuta* 
tien,  for  Ij  ui  tiquid  Xenon,  of  tedi  the  eqidBtarium  Ij  absoqMion  voctnun**  and  of  the  tiansiait 
aboorption  during  the  I2  photodiaaoGtation  foacdon  sequence,  whidi  we  wdl  <hacuas  in  mon  detail 
below. 


tUNELENGTH  (im) 


energy  (I0’em‘l 


ni.  4.  TfuHsieHt  ekctnmk  ^mctnm 
fiom  mokadar  tfynaimcsjbr  I2  photo- 
^Sssodaikm  bi  liquid  Xenon,  made  up 
of  tmnsUonsfrom  the  ground  suite  to 
die  exeded  A,  B,  and  B"  mtes.'''* 
The  congpaed  eguttihrtum  room  tem¬ 
perature  spectrum  Is  shorn  at  up  widi 
the  gat  phase  /j  band  conurur  as 
measured  by  T^t^udsen^*  dwm  as 
dots.  Briow  are  the  conputed  uan- 
slent  spectra  fiom  0  u>  900  ps  for  the 
reaction  sequence  beginning  aUh  pho- 
todtssodadve  excUation  to  the  A  state, 
then  either  /  aum  escape  or  sotvent 
caging  fidhwed  by  recombination  to 
form  a  highly  vAnokmady  excited  Ij 
mokctde  which  loses  ds  energy  u>  the 
sobrent.  ^ectra  are  sequeadady  plot¬ 
ted  with  6  X  lO'^i  vertical  rBs- 
pbcement. 


D.  X-ray 

Once  we  have  computed  the  molecular  dynamics,  i.e.  the  positions  of  atoms  as  a  function  of 
time,  we  can  then  compute  the  x-ray  diffraction  as  a  time  average  over  the  insuuntaneous  x-ray 
differential  scattering  cross  section.  The  basic  x-ray  scattering  equation  is 

in  which  is  the  differential  x-ray  scattering  cross  section  for  scattering  into  solid  angle  dfl 

dfl 

as  a  funclioo  of  the  photon  wave  (momentum)  dunge  vector  k  and  time  t;k*kf-k|isthe 
diffisienoe  between  the  wave  vector  kr  of  the  scattered  x-ray  photon  and  the  wave  vector  k^  of  the 


incident  idioton  in  ndiidi  Ikil  -  2ir/Xi.  Xi  being  the  incident  wevelength  end  the  dhection  of  kj 
being  the  incklmt  idioton  propegedon  direction;  e  is  the  diaiiB  on  die  etectron,  m  is  die  dectron 
mass;  c  dw  sp^  of  r  r«»eVsK^  '»  dw  risssirbl  adtas  of  die  electron  and 

t}  »  7.9398x10*^  cm^,  (l4€DS^)/2  is  the  factor  oonrecdia  for  the  uw  of  wqxiiMized  x-cap  pho* 
torn,  e  b^  the  scatterii«  angle  between  kf  and  ki;  nd  d^  p(r,t)  is  the  prob^Utjr  for  flni^  an 
dectron  in  volume  elenmit  d^  atpositioo  randthnet. 

Taking  into  account  the  molecular  synunetrp  and  the  adecden>7  of  motecutv  dignment  by  die 
pump  light,  tfid  using  standard  atomic  structure  factors,^  we  compute  as  illustrsfod  in  Fig.  S  the 
transient  x-ray  diffiraction*  from  the  I2  photodissociation  reaction  during  the  first  half  picosecond  of 
the  reaction  sequence.  Note  the  deviation  (which  is  small)  from  cyiindricd  symmetry  produced  by 
the  diotosdection  alignment.  On  the  experimentd  side,  in  collaboration  with  Mouiou  and  William¬ 
son  at  Rodiester,  we  have  recently  produced  x-ray  pulses  in  the  hundreds  of  picosecond  range  by  a 
light  to  electron  to  x-ray  conversion  process.^' 


diffractiim  for  1 2  dis- 
solved  In  (ydohexane, 
calculated  from 

molecular  cfynamics 
for  the  first  0.5  ps  of 
the  ti  photothssocia- 
tkm  reaction  se- 
tptence.^  The 

diffiuctkm  is  confut¬ 
ed  fitr  Titanium 
x-rays  (2.75  k).  The 


t^ffhtctkm  intensity, 
the  center  of  the  circle 
is  scattering  at  0*  and 
the  outer  edge  is 
scattering  at  90*. 


IV.  TRANSIENT  ELECTRONIC  ABSORPTION  SPECTROMETER 


In  order  to  check  our  computed  dynamics  against  nature's,  we  have  assembled  a  picosecond 
transient  dectronic  absorption  spectrometer,^^' **'*2  shown  in  Fig.  6.  A  mode-locked  Ar^  laser 
synchronously  pumps  two  dye  lasers,  producing  two  tunable  streams  of  picosecond  pulses.  Since  the 
energy  per  pulse  is  in  the  rU  range,  in  order  to  achieve  sufficient  excitation  our  sample  region  is  a 
smdl  confocd  volume  at  the  focus  of  microscope  objectives  in  both  the  pump  and  probe  beams, 
whidi  are  cdiinear  but  travdiing  in  opposite  directions.  The  reaction  solution  is  rapidly  pumped 
through  the  confocd  volume  inside  a  smdl  capillsry.  Because  the  cross  sections  for  pumping  and 
probing  are  both  severd  orden  of  magnitude  smdier  than  those  of  the  strong  absoiWs  used  for 
most  picosecond  experiments,  we  devdoped  a  triple  modulation  acheme‘^'“  using  inexpensive 
radio  amateur  gear  for  high  frequency  driving  and  detection.  This  modidation  system  uses  the 
disoovery  by  Heritage^^  and  Levine  and  Bethea^' ^  that  noise  in  Ai*  laser  systems,^' ^  in 
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V.  AN  EXAMPLE  REACTION:  h  PHOTODISSOaATION,  SOLVENT  CAGING,  RADI¬ 
CAL  RECOMHNATION,  AND  VIBRATIONAL  DECAY 

Even  tbough  little  is  yet  known,  it  fii  oar  btffief  Out  solution  phase  reaction  miriecular  dynam¬ 
ics  can  and  will  be  understood.  The  propw  apiKOKli  reqfuires  a  comMnation  of  theoretical  and  exper- 
imentM  Mduuques.  To  this  end  we  have  Mready  developed:  i)  the  theoretical  and  computational 
toois^  to  calculate  the  molecular  dynamics  of  Mlution  reactions;^'  ^  ii)  the  tools  to  calculate  from 
molecular  dynamics  the  interface  to  the  exp^ment,  i.e.  transient  spectra,^*  and  transient  x-ray 
diffraction;*  and  iii)  the  use  of  the  experimental  tool  of  picosecond  transient  electronic  absorption 
spectroscopy,^^  as  described  below. 

With  these  tools  we  have  studied  the  sample  case  of  the  reaction  sequence  in  various  solvents 
of  I2  photodissociation,  solvrat  caging,  radical  recombination,  and  vibrational  decay.  The  I2  reaction 
is  a  case  of  a  diatomic  involved  in  a  simple  reaction  for  whidi  the  electronic  spectra  of  the  reactant 
and  that  of  the  product  (the  same  molecule  in  this  special  case)  are  sufficiently  undNStood  that  tran¬ 
sient  electronic  spectra  can  be  used  to  decode  the  dynamics. 

A.  Theoretical  Moleailar  Djmmks 

The  first  detoministic  theoretical  study  of  the  molecular  (hmamics  of  solution  reactions  was  by 
Bunker  and  Jacobson,^  who  oomputed  the  dassical  tiqjectories  for  I2  in  OCU  solvent  represented 
26  spherical,  structurdeas  particles  in  a  specular  cube.  MurreU,  Stare  and  Dsmmel^  modeOed  the 
phototfissociatioa  of  I2  in  dense  inert  geses,  I2  plus  22  gm  atoms  in  a  spherkd,  soft-wided  oontainre. 

IMng  the  codes  we  have  developed^  and  the  array  processor,!’^  we  have  oompuied  the  molec¬ 
ular  dynamics  for  the  I2  reaction,  begtaudng  whh  phaiodiaBociadoa  through  the  A  state,  in  aevred 
dUBstent  solvents,  indwHog  ethyleoe  Rbwol,  gpclehanaae,  cartwn  tamcMorkle,  and  liquid 
xenon.  We  see  dramatic  (fiffisrences  in  the  molecalar  dynaiaics  among  Rm  dffbrent  solvents. 

The  oonduskm  of  all  three  molecular  dynamic  stucfies,  Bunkre  and  Jacobson,  Mmell,  Stare, 
and  Dammei,  and  our  own,  is  that  direct  gem^e  recombination  is  usually  a  very  fast  process,  ovm 
within  a  few  picoseconds.  An  important  cavva,  and  a  weakness  in  these  theore^  studies,  is  that 
the  process  elwrcby  the  I  ato.  disoodad^  m  an  excited  stme  potential  surface  re^  the  ground 


state  aunaee  an  wMcii  tiite  taoomMna  is  not  imdtestood,  and  b  tiierrfore  toidtBd  in  tease  cal- 
cidatkxK  Iqr  arMtiary  assumptions  which  may  be  incorrect.  If  so,  tee  red  time  for  geminate  recom¬ 
bination  may  be  longer  than  the  few  incoseconds  calculated.  These  calculations  suggest,  as  do  the 
non-molecular  dynamic  theoretical  computations  of  Nesbitt  and  Hynes, that  the  time  scale  of 
the  I]  reaction  sequmce  may  be  dominated  by  the  vibrational  energy  decay  step. 

B,  Theontkal  Transient  ^ectra 

From  the  computed  molecular  dynamics,  using  the  techniques  outlined  above,  we  calculate  the 
transient  electronic  absorption  spectra,  shown  in  Figs.  4  and  7. 


ETHYLENE  GLYCOL 


CYCLOHEXANE 


CARION  TETHACNLORIDE 


Fig.  7.  Theoretical  transient  electronic  ^tectra  for  I2  photodissodation  in  various  solvents  as  computed 
from  molecular  dynamics.^'^ 


As  can  be  seen,  the  differing  dynamics  among  the  solvents  are  reflected  in  differing  transient  spec¬ 
tra.  The  time  scale  of  the  dynamics  and  spectral  change  increases  strongly  in  the  order  ethylene 
glycol,  cyclohexane,  CG4,  i.e.  in  decreasing  order  of  strength  of  solvent  association.  Thus  the  two 
absorption  peaks  which  move  in  from  the  infrared  and  ultraviolet  evolve  back  toward  the  original 
single  visible  peak  much  more  sloaly  for  the  weakly  associated  solvents.  In  addition,  in  the  14,000 
to  17,000  cm~'  range  of  interest  here,  for  higher  energy  probe  photons  the  time  scale  of  the 
observed  transient  absorption  should  increase,  as  it  takes  longer  at  higher  probe  energy  for  the 
infrared  absorption  peak  to  arrive  and  to  pass  by. 
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C.  Ejqmimentat  TransieHt  Spectra 

The  earliest  picosecond  experimental  results  on  I2  were  by  the  Eisenthal  group,^i<32  who  meas¬ 
ured  the  tranaent  electronic  ab^ption  spectra  after  excitation  at  530  nm.  Decay  times  of  —70  ps 
for  I2  in  hexadecane  and  —140  ps  in  CQ4  were  observed.  Subsequently  these  studies  were 
extended  by  the  Eisenthal  group  to  I2  in  aromatic  solvents^^  which  are  telieved  to  form  complexes 
with  I  atoms,  and  by  Langhoff^^  who  observed  I2  photodissociation  in  several  weakly  associated 
liquids,  finding  decay  times  in  the  —100-150  ps  range.  More  recently,  Kelley  and  Rentzepis^^  have 
observed  I2  photodissociation  in  fluid  and  liquid  Xe  with  a  decay  time  of  —40  ps,  as  well  as  in  CCI4, 
and  similar  experiments  have  been  carried  out  by  the  Peters  group.^^  All  of  the  above  studies  used 
second  harmonic  Nd  pump  light  at  —530  nm,  which,  at  least  in  the  gas  phase,  results  in  excitation 
largely  to  the  bound  B  O^OW  state  which  presumably  predissociates,  but  may  absorb  another  pho¬ 
ton  in  the  meantime. A  delay  of  —20  ps  between  excitation  and  the  minimum  in  the  absorption 
curve  was  attributed  to  absorption  from  the  B  state.^*'^**^’  Thus,  the  linking  of  these  experimental 
transient  spectra  to  the  molecular  dyrutmics  of  the  I2  photodissociation  and  recombination  reaction  is 
made  difficult  by  the  probable  presence  of  the  additional  processes  of  B  state  predissociation  and 
absorption.  We  have  used  a  different  experimental  approach,  as  described  above,  to  pump  into  the 
dissociative  A  state,  thus  avoiding  the  problem  of  predissociation.  An  experimental  caveat  which 
should  be  remembered  for  all  such  Iodine  experiments,  ours  included,  is  the  well-known  propensity 
of  both  1  atoms  and  I2  molecules  to  form  complexes,  which  could  be  involved  in  the  observed  tran¬ 
sient  spectra.  Fig.  8  shows  our  experimental  measurements  of  the  transient  electronic  absorption 
during  the  I2  photodissociation  reaction  in  different  solvents  at  two  different  probe  wavelengths. 
When  the  theoretical  (Fig.  7)  and  experimental  (Fig.  8)  spectra  are  compared  as  a  function  of  time, 
of  solvent,  and  of  probe  wavelength,  there  is  gratifyingly  good  general  agreement. First,  as  our 
molecular  dynamics  theory  predicts,  the  time  scale  for  the  evolution  of  the  absorption  increases 
across  the  sequence  ethylene  glycol,  cyclohexane,  CQ4.  Second,  the  time  scales  are  at  least  semi- 
quantitatively  similar  in  theory  and  experiment.  Tliird,  the  observed  spectra  for  a  particular  solvent 
do  increase  in  time  scale  as  the  probe  photon  energy  is  increased,  again  as  predicted  by  theory.  This 
gives  hope  that  for  the  first  time  a  mqjor  portion  of  the  molecular  dynamics  of  a  solution  reaction 
may  indeed  be  understood. 

D.  Interpretation  of  Molecular  Dynamics 

It  would  not  be  truly  satisfying  to  only  describe  the  molecular  dynamics  of  solution  reactions  in 
terms  of  computer  calculations,  no  matter  how  well  they  agree  with  experiment.  We  would  like  a 
simple  physical  model  to  unify  the  observed  results.  In  this  light  we  offer  the  following  hypothesis 
which  fits  the  presently  known  theoretical  and  experimental  data.'*)  Let  us  consider  the  solvent  as  a 
bath  which  is  driven  by  the  motions  which  the  reactant  molecules  must  go  through  to  become  pro¬ 
ducts.  If  we  think  of  the  reaction  in  the  frequency  domain,  we  must  drive  the  bath  with  the  fre¬ 
quencies  of  the  atomic  motions  by  which  the  reaction  occurs.  To  a  first  approximation,  we  can  take 
a  linear  response  view  and  just  consider  the  natural  frequency  spectrum  of  the  unperturbed  sdvent 
as  the  frequence  response  of  the  solvent  to  being  forc^  by  the  reacting  molecules.  We  take  this 
spectrum  as  approximately  the  pure  solvent  velocity  spectrum,  i.e.  the  power  spectrum  of  the  veloci¬ 
ties  (or  equivalently  the  Fourier  transform  of  the  velocity  autocorrelation  function)  which  for  a  har¬ 
monic  system  would  be  the  normal  mode  density.^  For  the  solvents  of  interest,  these  spectra  are 
shown  in  Fig.  9.  This  simple  picture  is  related  to  the  more  sophisticated  and  fully  developed  treat¬ 
ment  of  solution  reaction  molecular  dynamics  recently  proposed  by  Brooks  and  Adelman.^i*^^ 

The  frequency  range  of  interest  for  the  bulk  of  the  I2  reaction  time  history  is  that  of  the  I2 
vibrational  motion  as  it  decays  to  the  solvent,  i.e.  the  1(X)  to  2(X)  cm~'  range.  As  Fig.  9  shows, 
strongly  associated  liquids  like  ethylene  glycol  have  a  high  spectral  density  in  this  region  because  of 
hindered  translational  and  rotational  motion,  and  we  would  thus  expect  the  reaction  sequence  to 
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proceed  swiftly  in  such  a  solvent,  as  seen  in  both  our  molecular  dynamics  caloilations  reflected  in 
Fig.  7  and  our  experiments  as  shown  in  Fig.  8.  Cyclohexane  has  a  lower  spectral  denaty  and  we 
would  expect  a  slower  evolution,  as  observed.  Finally,  CCU  has  the  least  spei^  density  in  this  fre¬ 
quency  range,  and  we  would  expect  it  to  be  even  slower,  again  as  observed  theoretically  and  experi¬ 
mentally.  Thus  velocity  spectrd  density  may  be  a  unifying  theme  in  the  imdeistanding  6[  these 
reactions  in  solution. 


Fig.  8.  E)q)erimental  tran¬ 
sient  electronic  absorption 
spectra"^'*  for  Ii  photodissoci¬ 
ation  in  various  solvents  at 
two  differetu  probe 
wavelengths.  Note  that  the 
time  scale  of  the  observed  ab¬ 
sorption  change  decreases  with 
solvent  association  strength 
and  increases  with  probe  pho¬ 
ton  energy  as  expected  from 
the  theoretical  spectra  shown 
in  Fig.  7. 


Fig.  9.  Velocity  spec¬ 
tra  for  various  sol¬ 
vents  as  computed 
from  nudecular 
dynamics.^^ 


VI.  DISCUSSION  AND  CONCLUSION 


One  of  the  most  important  unsolved  (fuestions  in  chemistry  .is  how  chemical  reactioas  take 
place  in  solution  in  toms  of  the  underlying  motions  oi  the  atoms.  We  have  already  in  place  basic 
theoretical  tools  needed  to  solve  this  question  for  a  wide  variety  of  inorganic,  organic,  and  biodiemi- 
cal  reactions.  We  can  compute  the  atomic  motions,  i.e.  the  molecular  dynamics,  for  reactant>soivent 
systems  invcdving  hundreds  to  thousands  of  atoms  for  times  on  the  femtosecond  to  nanosecond 
s^.2  We  can  calculate  from  the  dynamics  the  transient  infrared, Raman,^  and  decuonk^* 
spectra  and  x-ray  diffraction^  which  reflect  the  underlying  reaction  molecular  dynamics.  We  have 
now  b^un  to  test  the  predictions  of  these  theoretkal  techniques  against  experimentally  measured 
transient  spectra.  The  agreement  between  theory  and  experiment  for  our  first  test  case,  the  Ij  pho¬ 
todissociation  reaction,  is  sufficiently  good  that  we  are  encouraged  to  believe  that  we  are  at  last  able 
to  discover  the  molecular  dyiuunics,  the  microscopic  process,  by  which  a  real  solution  takes  place. 
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